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ISKANDER A. TAIMANOV AND SERGEY P. TSAREV 

Abstract. In this paper we consider the Moutard transformation [9] which is a two- 
dimensional version of the well-known Darboux transformation. We give an algebraic 
interpretation of the Moutard transformation as a conjugation in an appropriate ring 
and the corresponding version of the algebro-geometric formalism for two-dimensional 
Schrodinger operators. An application to some problems of the spectral theory of two- 
dimensional Schrodinger operators and to the (2 + l)-dimensional Novikov— Veselov 
equation is sketched. 



Introduction 

The Moutard transformation plays a fundamental role in projective-difFerential ge- 
ometry of surfaces and was extensively investigated by Bianchi, Darboux, Demoulin, 
Guichard and others. 

In recent publications [THl UHl EZ] we gave an application of the Moutard transforma- 
tion to the explicit construction of two-dimensional Schrodinger operators 

H = -A + u = -{dl + dl) + u{x, y) 

with fast decaying smooth rational potentials such that their i2-kernels contain at 
least two-dimensional subspaces spanned by rational eigenfunctions as well as a (2 -t- 1)- 
dimensional extension of the Moutard transformation which was able to produce explicit 
rational blowing-up solutions to the Novikov- Veselov (NV) equation with fast decaying 
smooth rational Cauchy data. 

Purely algebraic constructions of the aforementioned papers imply that one should 
look for an algebraic interpretation of the Moutard transformation as a simple transfor- 
mation in the ring of partial differential operators. The one-dimensional case of Darboux 
transformation was studied in |7j. Here we give an analogue of their results; as it turned 
out, for the Moutard transformation one needs to recourse to a more complicated ring of 
rational pseudodifferential operators (Ore localization of the ring of partial differential 
operators). As a by-product, we obtain a more natural version of the Dubrovin-Krichever- 
Novikov formalism [S1|H] of algebro-geometric two-dimensional Schrodinger operators. In 
the last Section we briefly sketch our previous results [TBI IT7j . 
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First we give an account of the algebraic theory of the Darboux transformations. 
0.1. The Darboux transformation. Let 



be a one-dimensional Schrodinger operator and let lo satisfy the equation 

Huo = 0. 

The function uj determines a factorization of H: 

H = A^A, A=^^+v, A^ = ^+v, v=^. (1) 
ax ax 



UJ 



Indeed we have 

A = ( , , 

dx J \ dx J dx^ 



AT A I d \ f d \ d^ 2 

A' A^ { — + v \ = - — — + 



and the equation 



is equivalent to Huj = 0. If w is real- valued we have A* — A^ . 
The Darboux transformatioi^is the swapping of A^ and A: 



or in terms of 



u: 



H = A^A -^H = AA^, 



U — + —f u ^ 



It is easy to check by simple computations that 

Proposition 1. // ip satisfies the equation Hip = Eip with E = const then (p = Aip 
satisfies the equation Hp = Eip. 

Remark. In general the Darboux transformation is defined for any solution to the 
equation Huj — cuj with c = const (see, for instance, [7]). In this case it reduces to the 
transformation of H' = H — c for which H'uj = 0. 

0.2. The Moutard transformation. This transformation was invented by Th. Moutard 
for the hyperbolic equation of the form t/i^;!, — u{x, y)4' = in the context of local differ- 
ential geometry of surfaces and was extensively studied and used in many areas of the 
local surface theory [4] . Here we give an elliptic version [15j \W\ [17] of this transformation 
suitable for our purposes. 

Let H he a, two-dimensional potential Schrodinger operator and let w be a solution to 
the equation 

Huj = (-A -I- u)uj = 0. 
Then the Moutard transformation of H is defined as 

~ ujI + UJ?. 

H ^ -A + u-2Aloguj ^ -A-u + 2 „ ^ . 



^Sometimes it is also called the Crum transformation due to 3 . Burchnall and Chaundy called it 
transference |2]. 
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Proposition 2. If (p satisfies the equation Hip — 0, then the function 9 defined from the 
system 

M), = -c.2(^) , (co0)y ^ ('^) (2) 

\ljj y y \tU / X 

satisfies HO = 0. 

By this definition, if 9 satisfies ([2| then 

9 + t — const, (3) 

UJ 

satisfies ^ for any constant t. 

We shall use the following notation for the Moutard transformation: 

M^{u) = u = u-2Aloguj, M^{if) = {6*+ te C}. 

For one-dimensional potentials the Moutard transformation reduces to the Darboux 
transformation. Indeed, let u — u(x) depend on x only and uj = f{x)e'^'^. Then / 
satisfies the one-dimensional Schrodinger equation 

Hof-{-^ + u)f = cf 

and the Moutard transformation reduces to the Darboux transformation of Hq defined 
by/: 

H — Hq - —— > H — Hq — TT^ • 

oy' ay' 
li g — g{x) satisfies i/off = Eg, then Hip = with ip — e^yg{x). We derive from ^ 
that 9 = e^vh{x) satisfies H9 = Qiih^ y^TTl {jt ~ f ' i-^- ^"^^^ multiple of the 
Darboux transform of q: h = — ^-^-y=Aq where Hn—c = A is the factorization of Hn — 

c defined by /. The inverse Darboux transformation is given hy g — ^ ^ ^ (s + 'z^) ^■ 
Remark. We can rewrite ^ as 

(a-H^)^ = z(a-^)^, (9+^)^ = -(5-^)^ (4) 

which implies 

■ B ■ uj{9) = iuj ■ B ■ ijj^^{ip), u!^^ ■ d ■ oj{9) = —iuj ■ B ■ LO^'^{p), (5) 

where, as usual, z — x + iy, B ~ Bz — \{Bx — iBy), B = Bg = \{Bx + iBy). This 
representation will be used later. 

There is another two-dimensional generalization of the Darboux transformation called 
the Laplace transformation. Its relation to integrable systems was recently studied in 

m- 

1. The SPECTRAL CURVES OF ALGEBRAIC SCHRODINGER OPERATORS 

Since weakly algebraic two-dimensional Schrodinger operators are defined in terms of 
their spectral curves on the zero energy level we recall the definitions of the spectral 
curves of one- and two-dimensional Schrodinger operators. 
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1.1. One-dimensional Schrodinger operators. Let a one-dimensional Schrodinger 
operator 

commute with an ordinary differential operator L of order 2n + 1. Then the Burchnall- 
Chaundy theorem [2] guarantees that these two operators satisfy a polynomial equation 

Q(H, L) - 0, Q(A, E)^\^- P2n+iiE) = 0, 

where P2n+i{E) is a polynomial of degree 2n+l. This equation defines the spectral curve 
r of the operator H. Its points parameterize the joint "eigenfunctions" , i.e. solutions of 
the equations 

H^j = EtP, Li) = AV-. 

Moreover these functions are glued (after some normalization) into a function ip{P, x) 
which is meromorphic in P e F and after completing T by adding an infinity point 
E ^ oo the eigenfunction ip gets a singularity at oo: 

where = is a local parameter on T near the infinity point oo |10l [5]. Therewith 
it is said that the operator H is algebraic (or algehro-geometric). 

Example. The spectral curve of a constant potential. Let u{x) = c be a constant 
potential, then F is given by = i? — c, 

and the covering (A, E) ^ E ramifies at E = c where ^'±(2;, c) = 1. 

1.2. Two-dimensional Schrodinger operators. Given a partial differential operator 
H, one says that operators Li , . . . , L„ generate a commutative (mod H) algebra A if 
they satisfy the following commutation relations: 

[Li,Lj] = DijH, [Li,H] = DiH 

where Dij,Di,l < i,j < n, are partial differential operators. This definition was in- 
troduced in [5 in which they considered in detail the particular case when H is the 
two-dimensional Schrodinger operator (probably with an electromagnetic field) in two 
variables and introduced the following definition: 
• a two-dimensional operator 

H = dd + vB + u 

is called (weakly) algebraic if it is included in a nontrivial commutative (modiJ) 
algebra A generated by operators Li and L2 in two variables, the operators Li 
and L2 satisfy the polynomial relation 

Q(ii,L2) = 0(modi/) 

and to a generic point of the algebraic curve Q{Xi,X2) = there corresponds a 
fc-dimensional space of functions ^ which satisfy the equations 

Hijj = 0, Liijj = X.ip, i = 1,2. 

The dimension k is called the rank of A. 
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2. The conjugation representation 

2.1. The Darboux transformation. The Darboux transformation admits the well- 
known representation as a conjugation in the ring of pseudodifferential operators which 
implies interesting corollaries concerning algebraic operators. 

Let Hlj = and ^ = — ^ = ^ • Then in the ring of pseudodifferential operators in 
X we have 

ax ax 
where the functions uj and lo^^ are identified with the operator of multiplication by them. 
From that we conclude 

Proposition 3. 

H = A-H-A-^. (6) 

Fixing Lo let us denote by 

M = A-M- A-~^ 

the conjugation of a pseudodifferential operator M. It is clear that if L commutes with 
H , then L commutes with H: 

[H,L\^Q => [H,L]^0. 
However, given an arbitrary differential operator M, we have 

, d 

M -Lu^ Mbj + M' ■ — 

dx 

where Mlo is the function obtained by applying M to a; and M' is a differential operator. 
Hence 

M ^oj ■ — ■ Lu^^ ■ M -oj ■ { 

dx \dx 

UJ - ■ uj~^ ■ Muj • 1 I • + w • 4- • • M' ■ w"^ 

dx \dx J dx 

and_^e conclude that 

M is a differential operator if and only if • Afuj = A = const, which means 

Muj = Xuj. 



Let Huj — and let L be a differential operator of odd order which commutes with 
H. We assume that its order is minimal with respect to this property. Let, by the 
Burchnall-Chaundy theorem, 

g(i/,i) = 0, Q{E,X) = X^-P{E). 

Then we have (see [7]): 

(1) if Luj = Xuj and P{E) has no a multiple root at i? = 0, then H and L generate 
a commutative ring of differential operators and 



Q{H,L) = 0; 
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(2) if Lu ^ ■ Luj ^ const, then H and Lq — LH generate a commutative ring of 
differential operators and 

(3) if Luj = Xoj and P{E) has a multiple root at = 0, then the action of the 
Darboux transformation is inverse to one described in the previous statement: 
H and some operator Lq generate commutative ring of differential operators with 

Q-{H,Lo) = 0, Q-{E,X)^X^-E-^P{E). 

In particular, this implies that 

• The Darboux transformation preserves the class of algebro-geometric one-dimen- 
sional Schrodinger operators. Moreover it always preserves the normalization of 
the spectral curve. 

In [7j it is noted that some facts mentioned above were discovered by Burchnall, 
Chaundy and Drach. 

The first two cases are demonstrated by the following 

Example. The Darboux transformation of the constant potential. Let u — c ~ const 
and L = 

ax 

(1) Let uj = e^/^~°o^. Then v — — \Jc — cq, = 0, and the Darboux transfor- 
mation is even trivial. 

(2) Let u = i(e^/^^=" e"^/^^=") cos(^c - cqo;). Then 

Co - c 



"Vc — Co tan(Vc — cqx), 



^ cos2(^/c — CqX) ' 
2(c-co) 



cos^(v^c — CqX) 

If c — Co < 0, then u is not periodic. 

An example of the third case is easily derived from the following explicitly computable 
examples which are interesting in themselves. 

Example. Rational solitons via the Darboux transformation. Let u — and uj — x. 
Then 

1 1^2 

V = -, Vx = ^, U = 

X X'' a;^ 

The spectral curve T = {A^ = E} is transformed to F = {A^ = E^} and 

V i^/ExJ 

(we normalize it by condition tp « ^iVex ^ _^ The iterations of the Darboux 
transformation initially applied to the trivial potential u = give all rational solitons 
discovered in [T]. The spectral curve of the potential 

_ n{n+l) 
2x'^ 

obtained after n iterations, is given by the equation A^ = The spectral curves of 

these are singular: topologically they are spheres but at i? = they have the following 
singularities: 
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any rational function f on r„, the spectral curve of m„, which is holomorphic near 
E = satisfies the condition 

f = /"' = • • — /('"-^^ = 0. (7) 
This is easily explained by the normalization mapping C ^ F = {A^ = ^2n+i| ^JiicJi 
has the form 

{X^t'^''+\E ^t^). 

Indeed any rational function / which is holomorphic at _E = is a ratio of polynomials 
gl^'^j such that Q(0, 0) 7^ and any such function written in terms of t satisfies the 
conditions ([7|. 

2.2. The Ore localization. In order to give an analogous interpretation of the Moutard 
transformation as a conjugation we first recall some basic facts from the Ore theory of 
localization of noncommutative rings and introduce the rings F(dy) and F{dy)[dx] which 
we shall use. 

We recall that a ring R is an algebra with two operations, the addition and the 
multiplication, such that i? is a commutative group with respect to the addition and the 
distribute laws 

a{b + c) = ab + ac, (a + b)c — ac+ be for all a,b,c Cz R 
hold. It is said that a ring R is regular if it satisfies the Ore conditions: 

• for a ^ 0, 6 ^ there exist r 0, s ^ such that ar — bs; 

• if a6 = ac or ba — ca for a ^ then b = c (so it has no zero divisors) . 

Ore showed in [2] that 

any regular ring can be embedded as a subring into a non-commutative Reld (skew 
field of fractions ) 

as follows. Let us consider the set S of all formal fractions 



ab 



Note that due to non-commutativity of R one shall always keep in mind that in our 
construction the denominators b^^ are always on the right of the numerators a! It is 
easy to propose an analogous construction with denominators standing on the left, but 
the resulting skew field will be isomorphic to the field we are to construct. 
We say that two such fractions are equal: 

if and only if 

ar = cs 

where r and s satisfy the equality (see the first Ore condition) 

br = ds. (8) 

It is easy to show that the equality is independent on the choice of r and s satisfying the 
last equality. Then the addition is defined by 



/ a\ / c\ f ar + cs\ f ar -\- cs\ 

\b) ^ Kd) = [-1^) = \-^) 
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The multiplication is given by 

ib) id) ^ ijhi) 

bt = au. 
The unit element is defined as 

© = '• 

this definition is independent on the choice of a and with such a law of multiplication we 
have (f) (f ) = 1. 

We say that S (with these operations) is the Ore localization of R. One can routinely 
check all the usual properties of the addition and multiplication operations defined on 
S, which becomes a skew field. This means that we can use expressions like A^^ for 
arbitrary elements of this skew field. All operations in S given above are constructive, 
unlike the case of the ring of formal pseudodifFerential operators defined as infinite series 
customary in solitonics. 

For a commutative ring R without divisors of zero the Ore localization coincides with 
the standard localization. Let us introduce the main example with non-commutative 
ring R which we shall use. 

Example. Let F — k(xi, . . . ,Xn) be a field formed by all k- valued functions of the 
form 

P(^Xi , . . . , Xji^ 
Q{xi , . . . , Xji) 

where P and Q are analytical functions of . . . , a:„ or even formal power series in these 
variables. We consider two cases which are 

1) k = K and a^i, . . . , a:„ G K; 

2) k = C, n = 2m and (xi, . . . ,a;„) = {zi,zi,. . . ,Zm,Zm), 

and, for simplicity and our needs, we even assume that n = 2, xi — x, and X2 ~ y- Let 
us consider _F as an operator algebra R which acts on itself by multiplications 

/(5) = /5, 

take the partial derivative operators dx — > c^y = and consider the operator algebra 
F\dy\ which is generated by elements of R and by the operator dy that acts on F. The 
ring is noncommutative because 

[dyj]^dyf-f-dy = ^^fy for feF. (9) 

It is straightforward to check 

Proposition 4. F[dy] and F[dx,dy] satisfy the Ore conditions. 

Let us take the Ore localization F{dy) of F[dy]. We have 
Proposition 5. For any f E F we have 

W-i fl = -g-i . ^ . g-i 

l^y ^ J i ^y ^y ■ 



where r and s satisfy (jSj). 



where 
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Proof. Let us take the equality ([9|, multiply every its term by dy^ from both sides 
and obtain [/, dy] — dy^ ■ fy ' dy^ ■ Proposition is proved. 

Let us consider F{dy) [dx] , the ring of formal differential operators in x with coefficients 
from F{dy). This ring is embedded into its Ore localization F{dx,dy). 

We need to define the commutation of d^^ with elements of F{dy)[dx]- Since 

[d-\f] = ^d-'[dxJ]d-\ 

it is enough to define the commutators of / G F{dy)[dx] with dx- Here we do that even 
for the Ore localization ring F{dx, dy). 

If L e F[dx, dy] is a differential operator in x and y with coefficients from F, then its 
derivative in x (i.e. the operator with coefficients being the derivatives of the respective 
coefficients of L) may be defined by the commutation formula 

Lx = [dx,L]. 

Let us take this formula as a definition of the derivative of any element of P e F{dx, dy) 
w.r.t. X. Using the Leibnitz identity 

{LM)x = LxM + LMx, L, M e F{dx, dy), 
we derive the following proposition by straightforward computations. 
Proposition 6. Given P ^ M ■ L^^ with L,AI £ F[dx,dy], we have 

P. [dx,P] = MiL-% + MxL-' = -ML-'LxL-' + MxL-\ 

Corollary 1. If P £ F{dy), then Px = [SjjjP] G F{dy), i.e. Px contains no derivations 
in X. 

2.3. The Moutard transformation. Let us consider the general Moutard transforma- 
tion which is applied to an operator: 

H = drds — u(r, s) 

where 

r = X, s — y, x,y CzM. (the hyperbolic case), 

or 

r — z, s — z, z G C (the elliptic case). 

Via the formulas 

{LUe)r ^ -iU^ , {iU0)s^LO^('^) , (10) 

it relates solutions ip and 6 to the equations 

Hip = 0, He = o 

where 

H = drds — u, u = u — 2drdsloguj = —u + 2 — 

Let us consider the Ore localization F{ds) of the ring of differential operators in s and 
the ring F{ds)[dr] of differential operators in r with coefficients from (noncommutative) 
ring F(ds). Given lu, a solution to the equation Hlu — 0, we consider the differential 
operators A,B G P[9s]: 

A = uj^^ ■ ds ■ ijJ, B = u ■ ds ■ uj~^ 
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and their "ratio" , the operator fl of the form 



LO ■ ■ UJ 



= A-^B. 



We denote by M and M the following formal operators of the first order from F{ds)[dr\: 

M = d-^ ■H = dr-d-^ -u, M = ■ H = dr - ■ u. 

Theorem 1. Given lo satisfying Hcu ~ {drds + u)lo = and the Moutard transform u 
of u (we assume that the transformation is generated by lo), the operators M and M are 
conjugated in F{ds)[dr] by fl: 

Remark. It is easy to check that H and H are not conjugated by fi: H ^ ■ H ■ il^^. 
Proof. First we expose some auxiliary facts. 



Lemma 1. 



n = A^^B =1 97^ • LUs- 

LU 



(11) 



Proof of Lemma. It is enough to show that B = A[l — ^ ■ ■ Ug). Let us write 
down the right-hand side of this equality as 

2 

A — u}^^ ■ ds ■ ijj— ■ dZ'^ ■ uJs — ■ ds ■ uj — 2uj^^ujg = 

LU 

= uj^^{ujs + uj ■ ds) — 2uj~^uJs — —uj^^ujs + ds = B. 
This proves the lemma. 

Corollary 2. 



(A-'Blr =1 d;'-LU, 



2uJr „_i 2 

• -UJs ■ -LOr 



(12) 



Now we are ready to prove Theorem[l] It to enough to establish the equivalent equality 

H -n^ ds-n - d^^ ■ H, (13) 



I.e. 



{drd, -u) ■ A-^ ■ B ^ ds ■ A-^ ■ B ■ d'^ ■ H. 
By Lemma [l] the left-hand side of ( 14 1 is 

{drds -u)-A-^-B^ds \{A-'^B).r + A-'^B- d, 
= ds ■ 



(14) 



i---d:^-u,, 

LO 



2uJr 2 

■ ds -UJs ■ ds -LOr 

UJ'^ LU 



.■ll---d7^-, 



2uJr 



■d-'-i 



2 

LU 

2uu. 



2L0rL0s 



LO 



2 -ds-ds^ -UJs - 



^ds ■ A-^B ■ dr 
l-'^.d-'-u:. 

LO 



LOrs • dsd ^ ■ LOrs 

LO 



ds ■ A ^B ■ dr+U 0,-LUs 7, 1 TT- ■ O^-LUs 
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= ^ ■ Os -^s ^ ■ Og • H n ^ ] ■ Os -^rs 

\ bJ^ LO"^ UJ'^ \ LW^ J UJ 

, « /(-Id O , 2^ «-l 2uJrUJs , ^UJr^^s 



• (ucj) -2u + ds- A-^B ■dr + u = 

In the right-hand side of ( 14 1 we have 

d, ■ A-^B ■ ■ {drds - u) = 9, • A-^B ■ dr - d, ■ (^1 ~ ^ ■ d''^ ■ oj,^ d''^ 



u 



d, ■ A-^B -dr-u+i-] -87^ -LUs-d^^ ■u + (-] -ds- dZ^ ■ ujs ■ ■ u : 



= (9, • A-^B ■dr-u-'^ -d-^ -LOs-d-^ ■u+ — ■ d-^ ■ u. 

Let us apply Proposition [s] to the third term ■ LOg ■ = uj ■ d^^ — d^^ ■ lo in the last 
formula and finally derive 

■ A^^B ■ dr - u-"^ ■ iuj ■ d-^ ~ d^^ ■ uj) ■ u+ — ■ d-^ ■ u = 



dg ■ A ^B ■ dr — u + ■ ■ oj ■ u, 



which coincides with the left-hand side of ( 14 1 . 
Theorem is proved. 

Let us assume that is a (weakly) algebraic operator, i.e., there are differential 
operators Li and L2 such that 

[Li,L2] = DoH, [Li,H] = DiH, [L2,H] = D2H (15) 

and there is a polynomial Q in two variables with constant coefficients such that 

Q(Li,L2) = (modi/). 

By applying the Euclid algorithm (division with remainder in the ring of formal linear 
ordinary differential operators F{ds)[dr]) we obtain -Ri,i?2 G Fids) such that 

Li-Qi-M = RieF{ds), L2-Q2-M = R2eF{ds) (16) 

with Q, e F{ds)[dr]. 

Theorem 2. (1) [Ri, R2] ^ 0; 

(2) Q(i?i,i?2) -0; 

(3) [i?i,M] = [i?2,M] = 0. 

Proof. 



(1) Using (151 we compute that 

[i?i,i?2] - [Li - QiH,L2 - Q2H] ^Si-H, SiQ F{ds)[dr]. 

However in the left-hand side we have an element of F{ds)[dr] which contains no 
derivations in r which implies that 5*1 = 0. 
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(2) The equality Q{Li,L2) = (modi/) means that Q{Li,L2) = U ■ H in F[dr,ds] 
which, by (15 1, imphes 

Q{Ri,R2)^ S2-H eF{ds)[dr]. 

However in the left-hand side we have an element of F{ds)[dr] with no derivations 
in r which implies S2 — 0. 

(3) We have 

[Ri,H] = [Li-QiH,H] = S3-H, 
[Ri,M] = [Ri, {ds)-^H] = {dsr^]H + {d,y^[Ri,H] = S4 ■ H = S5 ■ M, 
where Si G F{ds)[dr] for i — 3,4,5. By Corollary [l] this implies 

Af] = [R,,dr - • u{r, s)] - + • u(r, s)] G F{d,), 

and again as above we see that 5*5 = 0. The proof of the equality [R2, M] = is 
completely the same. 
Theorem is proved. 

In view of Theorems[T]and|2]we guess that the Moutard transformation should preserve 
the class of weakly algebraic two-dimensional Schrodinger operators. 

3. Applications of the Moutard transformation 

In [13 nil HZ] we gave a simple examples of fast decaying rational potentials for 
the two-dimensional Schrodinger operator which has a degenerated L2-kernel. These 
examples are constructed by using the Moutard transformation as follows. 

Main construction. Let 

Ho = -A = -A + uo 

be an operator with a potential uo{x, y) and let loi{x^ y) and uj2{x, y) satisfy the equations 

HqUi = HqUJ2 = 0. 

We take the Moutard transformations Ai^^^ and M^^^ defined by uji and W2 and obtain 
the operators 

i?i = -A + ui, H2^-A + U2 
where ui = ill^^ (uq), M2 — Mi^^{uo). By the construction, we have 

i?lAf^,(c^2) = 0, iJ2Mc..(t^l) =0. 

Let us choose some function 

Oi e M^,{uj2) 

and put 

02 = -— ^1 e Af^,K). 

UJ2 

These functions define the Moutard transformations of Hi and H2 and we obtain the 
operators H12 and _ff2i with the potentials 

-"12 = A^ei(wi), U2i^Me^{u2)- 
The following key lemma is checked by straightforward computations which we omit. 

Lemma 2. (1) U12 = "21 = u; 
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(2) For ipi — -gj and ip2 — we have 

Hill = Hii2 = 

where H = —A + u. 

We note that in this construction we have a free scalar parameter t (see (|3|) for the 
choice of 9i € M^-^ {^2)- This parameter can be used in some cases to build a non-singular 
potential u and functions ipi and '02 • 

For example if we apply this construction to the situation when uq — and uJi and 
UJ2 are real-valued harmonic polynomials 

3 

LJi = x + 2{x^ -y'^) + xy, UJ2 ^ x + y + -(x^ - y^) + 5xy, 

then for some appropriate constant of integration in 61 we obtain 

5120(1 + 8x + 2y+ 17x^ + ITy^) 

u — 

(160 + 4a;2 + 4y2 + iQx^ + 4x^y + IGxy"^ + Ay^ + 17{x^ + y'^YY 

5120|l + (4-i)z|2 



(17) 



and 



01 

•02 



(160-f |zP|2-f (4-i)z|2)2 
x + 2a;2 + xy — 2y'^ 



160 -I- 4x2 -I- 4y2 + I6a;3 + ^^^y + iQxy^ + 4y3 + I7(a;2 + y2y 
2x + 2y + 3.t2 + lOxy ~ Sy^ 



(18) 



160 + 4a;2 + 4^,2 _|_ iq^3 _|_ 4^,2^ _|_ iQxy^ + 4y3 + I7{x^ + y'^Y 
(here we simplify ipi and "02 by multiplying by some constant). 

Theorem 3 f |15lll6j ). The potential u given by |j7P is smooth, rational, and decays like 
l/r^ for r 00 (here r = -^/x^ -|- y^). 

The functions ipi and ip2 given by (18) are smooth, rational, decay like 1 /r^ for r — > 00 
and span a two-dimensional space in the kernel of the operator L — —A -|- u : L2{M.'^) 
i2(M'). 

If one takes appropriate two harmonic polynomials of the third order loi and UJ2 (cf. 
[m [16]) then one can construct the potential u and the functions ipi and "02 which 
are smooth, rational, the potential decays like l/r^ for r — > 00, the functions ipi and "02 
decay like 1 /r^ for r —> 00 and span a two-dimensional space in the kernel of the operator 
L^-A + u: L2(M2) ^ L2(K^). 

Remark. We guess that for every > by applying this construction to other har- 
monic polynomials one can construct smooth rational potentials u and the eigenfunctions 
■01 and -02 decaying faster than 

In [ini [TT] we used a time-dependent extension of the Moutard transformation con- 
structing explicit solutions of the Novikov-Veselov equation [TTJ [T^] 

Ut = d^U + B^U + 3d{VU) + 3d{VU) = 0, 

BV = dU. ^^^^ 

Some of our solutions show a very special behavior: the initial data for t = are smooth 
decaying rational functions of x and y; nevertheless for t > t^ > the solutions to the 
NV equation (19) blow up (become singular). 
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In particular the following solution U{z, z,t) of the Novikov-Veselov equation can be 
obtained using this technology: 

Hi 

with 

Hi = -\2{2Atx'^ + \2tx + 2Aty^ + \2ty + x^ -'ix-^y + 2x^ -2x'^y'^ -4,x'^y 
-2x^y^ - 60x^ - 3xy^ - ixy^ - 30x + y^ + 2y^ - 60y^ - 30y) , 

H2 = (Sx'^ + 'ix^ + 6a;2y2 + 3y4 4^3 ^_ 39 _ i2tf. 

This solution decays as at infinity, it is nonsingular for < t < = || and have 
singularities for i > = y|. 

For t^Tif, as we can see on Figure 1, the solution U{x,y,t) oscillates with growing 
amplitudes in the neighborhoods of the points P = (—1,0) and Q = (0,-1) since the 
denominator H2 vanishes at these points for t = T^,. The numerator Hi has zeros or order 
3 at the points P and Q so their respective neighborhoods are subdivided by smooth lines 
into 6 sectors with different signs of the numerator. The complicated behavior of the 
potential in the neighborhood of one of the singular points for i = is shown on Figure 2. 




y 05 y -0 



Figure 1: The potential [/ for i = j|. Figure 2: The potential U for t=^ 

in the neighborhood of the point (—1,0). 
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